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ABSTRACT 

We explore the possibilities for detecting pulsars that have ceased to radiate in the 
radio band. We consider two mod els: the model with hindered pa rticle escape from the 
pulsar surface [first suggested by I Ruder man fc Sutherlandlll975j | and the model with 
free particle escape (|Aronslll981l ; lMestellll999l) . In the model with hindered particle 
escape, the number of particles that leave the pulsar magnetosphere is small and 
their radiation cannot be detected with currently available instruments. At the same 
time, for the free particle escape model, both the number of particles and the radiation 
intensity are high enough for such pulsars to be detectable with the presently available 
receivers such as GLAST and AGILE spacecrafts. It is also possible that extinct radio 
pulsars can be among the unidentified EGRET sources. 
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1 INTRODUCTION 

The number of radi o pulsars in the modern catalogues al- 
ready exceeds 1500 ^Manchester et al"1l2005l ) and it is still 
increasing. Most of the observed pulsars belong to young 
neutron stars with the characteristic ages of td ~ 10 Myr. 
Since the lifetimes of radio pulsars are much shorter than 
the age of the Universe, up to 10 9 older neutron stars 
that have already ceased to radia te in radio band can ex- 
ist in our Galaxy (|Lipunovlll987h . The possibility of their 
detection is usually associated with the thermal radiation 
due either to internal stores of neutron -star thermal energy 
l|Yakovlev. Levenfish fc Shibanoy ]^999) or to accretion from 
the interstellar gas (|Colpi et al.lll998f) . Either way the elec- 
trodynamic processes in the neutron-star magnetosphere are 
assumed to cease to play a crucial role. However, the ces- 
sation of secondary-plasma production in the region of the 
magnetic poles at spin periods P ~ 2 — 3 s by no means 
implies that the electrodynamic processes become unim- 
portant. Thus, for example, it is well known that the sec- 
ondary plasma cannot be produced in the region of closed 
magnetic field lines only a t periods P > P cr ~ 10 5 B^ 3 s 
(|lstomin fc Mosvagin|[l995T) . where B 12 = B/W 12 G. There- 
fore, at P < P cr some of the closed magnetic field lines still 
remain filled with secondary electron-positron plasma. As we 
see, the critical period is long enough (more than a month), 
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so the electrodynamic processes in extinct radio pulsars play 
a significant role for quite a long time. Accordingly, parti- 
cles can be effectively accelerated in the region of open mag- 
netic field lines where a strong longitudinal electric field is 
generated. This acceleration results in observable radiation 
associated with curvature losses. 

Clearly, the total intensity of this radiation is related 
to the number of charged particles that fall into the region 
of strong longitudinal electric fields. This, in turn, directly 
depends on the work function of particle escape from the 
neutron star surface. Thus, analyzing the radiation from ex- 
tinct radio pulsars, we can draw a definite conclusion about 
the particle-escape work function and, hence, obtain inde- 
pendent information on the structure of the particle accel- 
eration region near the magnetic poles of radio pulsars. 

In this paper we discuss two basic models: the model 
with hindered particle escape from the neutron star sur- 
face (|Ruderman fc S utherland 19 75h and the model with 
free particle escape ( Aronsl [l981: Mcstcl 1999). We show 
that in the scope of model with hindered particle escape, 
the number of particles leaving the pulsar magnetosphere is 
small. So, their radiation cannot be detected with currently 
available instruments. At the same time, for the model with 
free particle escape, the number of particles and, hence, the 
total energy release are large enough for such extinct radio 
pulsars to be detectable with GLAST or AGILE satellites. 
Also, within the scope of free particle escape model, extinct 
radio pulsars can possibly be among the unidentified sources 
from the 3rd EGRET Catalog. 
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2 MODEL WITH HINDERED PARTICLE 
ESCAPE FROM THE PULSAR SURFACE 

The simplest model we use to estimate the radiation from 
extinct radio pulsars is the model with hindered particle es- 
cape from the pulsar surface. The effects of general relativ- 
ity for this model do not lead to any significant corrections, 
because they do not qualitatively change the structure of 
the el ectrodynamic equations in the particle-acceleration re- 
gion (Bcskm] |l999l ). Since the particle-escape work function 
for a surface of a neutron star in the Ruderman-Sutherland 
model is large, particles only fill the equatorial regions with a 
strongly curved magnetic field, in which secondary particles 
can still be produced. As a result, the longitudinal electric 
field outside the plasma (in particular, on open magnetic 
field lines) can be strong enough for effective particle accel- 
eration. In this case, the polar cap size is determined by the 
last closed magnetic f ield line on which seconda ry plasma 
can still be produced (|lstomin fe Mosvaginlll995l ): 

/ P \ 3/8 

Here R is & radius of a neutron star. 

Thus, the observed radiation of extinct radio pulsars 
is determined solely by the primary particles radiation. At 
that, the same curvature radiation mechanism as for the nor- 
mal radio pulsars can be considered. Recall, that curvature 
radiation is similar to synchrotron radiation, but the curva- 
ture radius a magnetic field li ne R c is used instead of the Lar- 
mor radius tl ~ m e c 2 -y /eB f Zhe leznvakovl [79771 ) . However, 
in contrast to young neutron stars, the condition for the pro- 
duction of secondary electron-positron pairs in the polar re- 
gions of extinct radio pulsars is violated. Therefore, in case of 
extinct radio pulsars the radiation should be expected only 
in high-energy part of the electromagnetic spectrum. We 
stress that this point distinguishes our analysis from other 
calculations of gamma-ray emission of normal radio pulsars, 
for which the entire se condary-particle spectrum should be 
taken into account fsee lHarding. Muslimov fe Zhanel (|2002l ) 
and references therein). 

We estimate expected flux from radio-quiet pulsars as 
the product of the energy losses by a single particle and the 
number of particles leaving the pulsar magnetosphere: 



F = 



ASN 
And 2 ■ 



(2) 



Here, N = dn/dt is the number of particles that leave a 
pulsar magnetosphere, AS is the energy losses by each of 
these particles, and d is the distance between a pulsar and 
an observer. Since the radiation mechanism in this case is 
curvature rad iation, the particle e nergy losses can be defined 
as AS. Here l|Zheleznvakovlll977V 
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~dl 
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is the power of the curvature radiation, R c ~ R 2 /R±^ is the 
curvature radius of magnetic field lines, and R± ta 10 4 cm 
(Eq. [TJ is the distance between the last closed magnetic field 
line and the magnetic dipole axis. 

It is well known that a strong longitudinal electric field 
exists only near the neutron-star surface. At rather small 
distance from a neutron star surface the longitudinal elec- 
tric field drops quickly even if there is no plasma flowing 



away in this region (jMuslimov fe Tsvgan|[l990h . Thus, we 
should consider two possible scenarios: the particle emits all 
its kinetic energy near the neutron star surface or it leaves 
the region of strong longitudinal electric fields virtually not 
losing its energy. In the former case, the particle energy 
losses AS are determined by the "radiation time" T ra£ j, and 
AS = mc 2 j. In the latter case, the energy losses are defined 
by the "escape time" (the time that a particle needs for 
leaving the region of the strong longitudinal electric fields). 
Generally, the minimum time At = min(Tesc, i" rac j) can be 
taken as At. The escape time is resc = Rc/c ~ 10 -3 s. 
For the radiation time we can write the following expres- 
sion T rac j = mc 2 ^y/(dS/dt), where the Lorentz factor of the 
accelerated particles can be estimated as 



eER± e(4TvpGjR±)Ri 



4iveR ± 2 B 
Pm e c 3 



(4) 



Here, Bo is a magnetic field on a neutron star surface. As a 
result, the radiation time can be written as 



7 rad 



R 4 c 5 
n 3 ulr e R 8 , 



(5) 



where pej = encj = — Bf2/(2-7rc) = B/(Pc) is the 
Golreich- Julian density, r e = e /m e c is the classical elec- 
tron radius, and uib — eB/m e c is the cyclotron frequency. 
Here we want to note that the radiation and escape times are 
of the same order: T rac j ps 10 _4 -Bi2 s, resc = Rc/c « 10 -3 s. 
Thus, within the accuracy of our estimate we can use either 
of these times. Finally, given the particle Lorentz factor, the 
characteristic energy of the curvature photons can be esti- 
mated as 



E ph « 2.9 x 10 3 Bi2 _1/2 J 



-3/8 



MeV. 



(6) 



We see that the expected energies are within the hard 
gamma-ray spectrum. 

As we mentioned above, within the scope of hindered 
particle escape model the particle production in the open 
field lines region is not possible. Still primary particles can 
appear in the region of strong longitudinal electric field. Here 
are two possible ways for this region to have a number of pri- 
mary particles. Firstly, the primary particle production can 
be associated with diffuse Galactic gamma-ray radiation, 
which causes the single-photon conversion of photons into 
electron-positron pairs. In this case the number of particles 
being produced nea r the magnetic poles can be estim ated as 
2 x 10 5 particles/s (jShukre fc Radhakrishnanll 19821 ). Hence, 
the expected flux (Eq. [2} is 



F w 0,3 x 10 _28 di"o 2 MeV/c 



where dio = d/10 pc. 

Another source of primary particles is a diffusion from 
the region of closed magnetic field lines filled with secondary 
plasma. The number of partic les leaving the pulsar magn e- 
tosphere can be estimated as l|lstomin fc Mosvagini ri995) 

at ujb 

As a result, the number of particles that diffuse into the 
region of a strong longitudinal electric field can be estimated 
as 



N = 



4-iT 2 m e jR 3 



4- 10 i8 P 



is D -4 particles 



(9) 
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The expected flux from the object (Eq. [2} is 
F » lO" 13 ^ 2 ^" 4 MeV/cm 2 s. 



(10) 



Thus, providing a validity of the model with hindered 
particle escape, the radiation intensity of extinct radio pul- 
sars turns out to be so low that the radiation power is not 
enough for the signal of neutron star to be detectable even at 
P ~ 2 — 4 s. The main reason is that the number of primary 
particles that fall into the region of a strong longitudinal 
electric field is very small, i.e., because the particle escape 
from the neutron star surface is hindered. 



3 MODEL WITH FREE PARTICLE ESCAPE 
FROM THE PULSAR SURFACE 

In this section we consider the model with free particle es- 
cape from the neutron star surface. Within this model there 
is only a slight difference between the charge density of the 
particles flowing away from the polar cap and the Goldreich- 
Julian density. General relativit y effects become sign i ficant 
near the ne utron star surface 1 Muslimov fc Tsyganl [l990l ; 
Beskinlll990Tl Since no particles are produced, the model of 



Mestell[|l999h . in which no secondary plas ma is generat ed, is 



appropriate in this case. The other model |Aronslll98lh can- 
not be realized for extinct radio pulsars because it requires 
a reverse flo w of secondary particles (for more details, see 
iBeskinll 19991 ). Here a short explanation should be done. The 
Mestel's model presumes that none of the magnetic field 
lines near the pulsar surface are preferred because of the ef- 
fects of general relativity. Therefore th e particles cannot be 
accelerated in this region (|Beskinlll990h . In other words, the 
longitudinal electric field that arises from the mismatch be- 
tween the plasma charge density and the Goldreich charge 
density decelerates rather than accelerates plasma particles. 

Nevertheless, the absence of regular acceleration does 
not imply that the plasma cannot fill the polar regions. Only 
effective particle acceleration is impossible. On the other 
hand, even at small distances from the pulsar surface r > ro, 
where 



ro 



1,8-10° 



M 

Mo 



(11) 



the general relativity effects become negligible. At such dis- 
tances the difference between the plasma charge density and 
the Goldreich charge density results in particle acceleration, 
at least on the half of the polar cap where the magnetic 
field lines deviate from the spin axis (i.e., those magnetic 
field lines for which the angle between the magnetic field 
line and the spin axis increases with distance from the neu- 
tron star). Below, we assume that particle acceleration be- 
gins only after the radius r — ro is reached. At smaller radii 
plasma on open magnetic field lines rotates rigidly with a 
neutron star. Finally, recall that in Mestel's model, as in 
Arons's model, the number of particles leaving the magne- 
tosphere of the neutron star is close to the Goldreich- Julian 
density near the polar-cap surface, i.e. 



N = TvR n G jc : 



10 30 fl 



ft 



particles/s. 



(12) 



Here, Ro = (QR/c)R is the polar cap radius. 



To define the electric potential, we write the Po isson 
equation in a rotating coordinate system (|Mestellll999T ) 



Aip = 4ir(pGJ - Pe) 



(13) 



Here, the Goldreich- Julian charge density and the particle 
charge density can be written as 



Pgj = 



Q.B _ Q.B 

2nc 2nc 

QB 

o e = c 

2ttc 



COS Ub, 



g(0) 



(14) 
(15) 



where 6b is the current angle between the magnetic field line 
and the spin axis, and #£ ' is the angle between the magnetic 
field line and the spin axis at the bottom of the acceleration 
region r — ro- 

To simplify our calculations, we consider axisymmetric 
case of an aligned rotator, for which all magnetic field lines 
are preferred. Using the expressions for the plasma charge 
density and the Goldreich charge density, we can write the 
following equation 



A(fi : 



9 nB(r ) 



(16) 



The boundary conditions can be defined as 



<p(r = ro) = 0; ip (r = ro) =0; ip (r = Rl sin 2 9) = 0. 

Here, -Rl = c/f2 is the radius of the light cylinder. As a 
result, for a zero approximation we have 

VBoR 3 

P-—-5-, (18) 



Rl' 



Thus. 



)Tf) 



(19) 



In the above expressions we used the fact that magnetic 
flux should be constant ("frozen" magnetic field lines), which 
yields B(ro)ro 3 = BoR 3 . The radiation and escape times can 
be estimated as i~ rac j and resc, the same way as for the model 
with hindered particle escape from the neutron star surface. 
Within the model with free particle escape T" rac j ~ i~esc at 

P = P break> where 



break 



2ttR 1 2r e \ 1/7 



iA 3/7 



0.95 s. 



(20) 



As we see, the break takes place at spin periods of greatest 
interest - when a neutron star just crossed the death line. 

The estimated energy of the curvature gamma-ray pho- 
tons leaving the magnetosphere of the neutron star E p h — 
hey 3 /R c « 8 ■ io^P" 13 / 2 erg « 0,05P" 13/2 GeV is within 
the gamma-ray range. Finally, the radiation intensity can be 
estimated as 



W : 



W 



5.8 • 10 



3i B\ 2 erg 



if P < P break 

2e 2 c 4 - (VLR\ 11 /ujbR\ 5 2m e c 3 

if P > P break . 



! )"(^y 



3R 



4 • 10 c 



Bl 2 erj 



The flux from such neutron stars near the Earth can be 
defined by the following expression 
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F : 



Wi erg 
laid 2 cm 2 s ' 



(23) 



where d is the distance to a neutron star, and ao = 
Rq/R is defined by the polar cap size Ro. Thus, for an 
oblique rotator the averaged flux from an extinct radio pul- 
sar is 

2Ro F photons 



Fnor 



2n R sin x E p h cm 2 s 



(24) 



Here \ is the angle between the magnetic and spin axes. 
As a result, the expected radiation fluxes from extinct radio 
pulsars within the scope of free particle escape model are 



Fi 



■ ■ w-*bU 13 r1 2/1 % c 2 e;^ 3 

p cm z s 
P < f'kreakj 

8 • IQ^BTRf^d^E'J 13 



p > p, 



break • 



(25) 



Since the expected fluxes are large, we conclude that for 
the distances d ~ 100 — 200 pc characteristic for the nearest 
neutron stars, such fluxes can most likely be detected with 
GLAST and AGILE satellites. 

On the other hand, since the radiation intensity strongly 
depends on the spin period of the neutron star, radio-quiet 
pulsars with long spin periods (i.e., those that have long 
ceased to radiate in the radio band or those with compara- 
tively long spin periods at the stage of radio emission) can 
no longer be detected. So there are three conditions under 
which extinct radio pulsars are detectable: the validity of 
the model with free particle escape from the neutron star 
surface, relatively small distances to such stars, and short 
spin periods of these stars. 

Figure [T] compares the threshold sensitivities of the 
modern gamma-ray missions (particularly, GLAST, AGILE, 
and EGRET) and the photon count rates for the model with 
free particle escape from the pulsar surface. We assumed the 
distance from a neutron star to the observer is 100 pc, ra- 
dius a star is 10 km, and magnetic field is 10 12 G. As we 
see, the closest extinct pulsars (if any of them exist in such 
a vicinity) can be detectable by the future missions such as 
GLAST and AGILE. 



4 CONCLUSIONS 

Thus, we have shown that extinct radio pulsars can be ob- 
served for some time as intense gamma-ray sources. For 
nearby neutron stars, they can be detected with the receivers 
installed on the modern gamma-missions. 

On the other hand, since the dependence on the spin 
period P of neutron star in the expressions for the radiation 
intensity is very strong, radio-quiet pulsars with sufficiently 
long spin periods (i.e., those which have long ceased to ex- 
ist as radio pulsars or those for which the spin period was 
relatively long even at the radio emission stage) still cannot 
be detected. Therefore, we can formulate three conditions 
for the possibility of detecting gamma-ray emission from ex- 
tinct radio pulsars: a validity of the model with free particle 
escape, a relatively small distance to such stars, and short 
spin periods. Providing all three conditions are satisfied, the 



detection of gamma-ray emission with the parameters de- 
scribed above would be a direct evidence for the free particle 
escape from the neutron star surface. 
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Figure 1 . Comparison of the sensitivities of present and future detectors in the gamma-ray astrophysics {Morsclli 20(3) and the predicted 
radiation intensities of extinct radio pulsars. 



